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Abstract: Thermostats models in space dimension d = 
1,2,3 for nonequilibrium statistical mechanics are con- 
sidered and it is shown that, in the thermodynamic limit, 
the evolutions admit infinitely many constants of motion: 
namely the intensive observables. 



I. THERMOSTATS 

Systems in nonequilibrium statistical mechanics have 
to be thermostatted, and the most realistic thermostats 
are infinite (i.e., very large) systems initially in equilib- 
rium. The present paper discusses a small (classical) 
system interacting with fairly realistic infinite (classical) 
thermostats of dimension d = 1,2, or 3. This can be 
attacked via recent difficult results on the time evolu- 
tion of infinite systems in dimension d = 1, 2, or 3. It is 
assumed that the initial equilibrium states of the ther- 
mostats are away from phase transitions. Some technical 
assumptions on the interactions are also made. The re- 
sult obtained here may then be expressed physically as 
follows: at any finite time each thermostat remains close 
to equilibrium in the sense that its global temperature 
remains the same, and this is also true for other inten- 
sive thermodynamic variables. If an infinite time limit 
were taken the situation would probably be quite differ- 
ent (and nontrivial only in dimension d = 3) but this is 
a hard problem, and not tackled in the present paper. 

The class of models that we shall investigate is when 
particles of a test system, in a container fio, and v other 
particles systems, in containers ill, . . . , interact and 
define a model of a system in interaction with v ther- 
mostats, if the particles in fii, . . . , f2„ can be considered 
at fixed temperatures Ti, . . . ,T„. 

A representation of the system is in Fig.l: 




H 



x = (Xo, Xo, Xi, Xi, . . . , X„, X„) 
Fig.l: If d = 1,2 the 1 + v finite boxes flj n A, j = 
0, . . . , v , are marked Co, Ci, . . . ,C V in the first figure and con- 
tain No, Ni, . . . , N v particles, out of the infinitely many par- 
ticles, with positions and velocities denoted Xo, Xi, . . . , X„, 
and Xo,Xi,..., X„, respectively, contained in > 0. 

The second figure illustrates the special geometry that will be 
considered for d = 1,2, 3: here two thermostats, symbolized by 
the shaded regions, , 0,2 occupy half-spaces adjacent to Qq • 



From the point of view of Physics the temperatures 
in the thermostats are fixed. A natural model, often in- 
voked in the applications, [1], is to imagine the contain- 
ers Oj ;, j = 1, . . . , v, as infinite and occupied by particles 
initially in a Gibbs distribution with given temperatures 
and densities Ti, S±, . . . , T„, S u . 

To implement the physical requirement that the ther- 
mostats have well defined temperatures and densities the 
initial data will be imagined to be randomly chosen with 
a suitable Gibbs distribution 



Initial data: The probability distribution fiQ for the ran- 
dom choice of initial data will be, if dx d = II j=n ^ah*" 3 > 

the limit as A — > oo of the distributions on the configura- 
tions x G %(A) with Xj G A (see Fig.l), 



/j Q -^(dx) = const e dx 



ELo/W(*i) 



with H (x) 

a d ^S 1 

Pj ~ k B Tj ^ n l 

0, Aq G M will also be fixed 



(1.1) 
Uj(X-j)) and 



> 0, Xj G R,j > 0; the values (3q 



> 



The values /?o, An. bear no particular physical meaning 
because the test system is kept finite. Here A = (An, 
Ai, . . . A^) and T = (To, Ti, . . . , T v ) are fixed chemical 
potentials and temperatures, and A is a ball centered at 
the origin and of radius ro . The Kj (Xj ) , Uj (Xj ) are 
kinetic and potential energies of the particles in fij (see 
below for the conditions on the potentials). 

The distribution fiQ is interpreted as a Gibbs distri- 
bution /in obtained by taking the "termodynamic limit" 
A — > oo. At time we switch on the interaction be- 
tween the particles in fio and those in the thermostats 
£!j, j > 0. The measure /xo is n °t time invariant un- 
der the corresponding dynamics (existence of dynamics 
in infinite systems is not at all a trivial issue, as already 
revealed by the theory of the evolution in the infinite 
space [2] and as it will be discussed later) and we need 
to: 

(i) define the temperatures of the thermostats (which are 
outside equilibrium); 

(ii) prove that the "macroscopic" property of the ther- 
mostats of having given densities and temperatures re- 
mains when the system evolves in time. 

\def i 



If P j(/3,A;Ar 



logZj(/3, A) with 



OO „ 

N=0 



dx tv 



-P(-\N+K j {x N )+U j {x N )) 



(1.2) 



where the integration is over positions and momenta of 
the N particles in A n ilj then we shall say that (at least 
at time 0) the thermostats have pressures pj (/3j , Aj), den- 
sities Sj , temperatures Tj , energy densities ej , and poten- 
tial energy densities Uj, for j > 0, given by equilibrium 
themodynamics, i.e.: 

Pj(P,X) = J im Pj {Pj > A j > A ) 

A— >oo 



2 



k B T 3 



(1.3) 



j ~ dp, j j ' 



which are the relations linking density 5j, temperature 
= energy density e 3 and potential energy 

density u,- in a grand canonical ensemble and in absence 
of phase transitions in correspondence of the parameters 
(/3j,Aj), for j >0. 

Remark: (1) notice that the limit defining p., does not 
depend on the shape of flj and coincides with the usual 
definition of pressure in the thermodynamic limit in the 
sense of Van Hove, [3]. 

(2) As usual in Physics we could define density, energy 
density and temperatures in single configurations x as 



lim ( N i^(x) UjaJx) K jA Jx) \ 
n^oo V | A„ PI ilj I ' I A„ n Qj | ' iVj,A n (x)J 

provided the limit exists. 

(3) By the Birkhoff theorem applied to systems in the 
full space R d , the limits exist with probability 1 for 
any translational invariant infinite-volume Gibbs mea- 
sure (i.e. a DLR distribution, [4]). Moreover under an 
additional assumption of "extremality" the limits are al- 
most surely the same for all x. By suitable assumptions 
on the parameters f3j and Xj, stated later in this sec- 
tion, we shall see that the limits in Eq.(1.4) exist with fi 
probability 1 and are equal to the values in Eq.1.3. 

Time independence of the intensive observables (in 
particular those in Eq.(1.4)) is the central issue in this 
paper. Even if the evolution is defined with only Ho, 
i.e. no interaction between Qq and the thermostats so 
that no is time- invariant, yet, in general, one can only 
conclude that along "typical trajectories" the intensive 
observables are constant at countably many times (for 
instance at all rational times) . 

However under our assumptions on (3j and Xj (essen- 
tially absence of phase transitions) and in the interest- 
ing case when the interaction between fi and the ther- 
mostats is switched on then, by choosing the initial con- 
figurations with fi probability 1, we shall prove that the 
intensive observables keep the same initial value at all 
finite times. This justifies our terminology to call ther- 
mostat the systems Qj, j > 0. 

Hypotheses: In the geometries of Fig. 1 suppose: 

(1) Ho satisfies the DLR equations and that 

(2) the thermostats pressures pj(f3,X) are differentiable 
in (3, X at [3j,Xj, j = 1, . . . , v. 

It is essential that the "macroscopic" property of the 
thermostats, of having given densities and temperatures, 
remains when the system evolves in time. 

Evolution is defined via equations of motion: since 
we are dealing with infinitely many particles it will be 
defined by first considering the motion of the particles 



initially contained in some ball A keeping the particles 
outside A fixed. Such motion x — > S^x is called A- 
regularized: then we shall consider the limit as A — > oo. 

The regularization boxes A will be (for simplicity) balls 
A„ centered at the origin O and with radius 2"r ¥ ,, with r v 
equal to the range of the interparticle potential, and par- 
ticles will be reflected at the boundary of A„. The limit 
motion reached as n — > oo will define the thermodynamic 
limit motion. 

The A„-regularized equations of motion will be 



mX„, = -diU (Xo) -J2diU 0}j (X , X.O + ^Xq) 

j>0 

(see Fig.l) where: 

(1) the first label, j = or j = 1, . . . , v, refers (respec- 
tively) to the test system or to a thermostat, while the 
second indicates the derivatives with respect to the co- 
ordinates of the points located in the corresponding con- 
tainer and in the regularization box A„ (hence the labels 
i in the subscripts (j,i) have Nj d values). 

(2) The forces <P(Xo) are, positional, nonconservative, 
smooth stirring forces, possibly vanishing; the other 
forces are conservative and generated by a pair poten- 
tial ip, with range r v , which couples all pairs in the same 
containers and all pairs of particles one of which is lo- 
cated in Oo and the other in Clj (i.e. there is no direct 
interaction between the different thermostats). 

(3) Furthermore particles are repelled by the boundaries 
of the containers by a conservative force of potential en- 
ergy ip, diverging with the distance r to the walls as r~ Q , 
for some a > 0, and of range C r,,. The poten- 
tial energies will be Uj(Xj), j > 0, and Uoj(X-o, Xj), 
respectively denoting the internal energies of the various 
systems and the potential energy of interaction between 
the test system and the thermostats: 

^•(x)= E ^-?') 

qeXj (q,q')eX 3 ,qeA 

C/o, J (X ,X J )= E Vb-J) (!-6) 

The potentials ip, ip have been chosen j-independent for 
simplicity. 

(4) The equations are formally defined also in the 
phase space H of the locally finite configurations x — 
(...,%, qu... )gi 

x = (X ,X , Xi, Xi, . . . , X„,X n ) = (X, X) (1.7) 

with Xj C Oj, hence X C tt = U^ =0 Oj, and & e M d ; in 
every ball E(r') of radius r' and center at the origin O, 
fall a finite number of points of X. 

Infinite systems are idealizations not uncommon in sta- 
tistical mechanics. But we take it for granted that they 
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must be considered as limiting cases of large yet finite 
systems. This leads to several difficulties: one is im- 
mediately manifest if one remarks that the equations of 
motion Eq.1.5 do not even admit an obvious solution in 
H. 

Dynamics is well defined with ^-probability 1 be- 
cause if d = 1,2,3 the A„-regularized equations with 

def 

data x admit, with /zo-probability 1, a limit SfX = 
liniA^oo 5j A "'x for all t > 0: a precise statement is in 
theorem 4 below (proved in [5, theorems 6,7], for d = 1,2, 
and in [6, Theorem 1] for d = 1, 2,3). 

Since the Eq.1.5 are Newton's equations we shall call 
the model a frictionless thermostats model: this is to con- 
trast it with other thermostats models in which artificial 
"frictional" forces are introduced to make it possible for 
the system to reach a stationary state. In models with 
friction entropy production (generated in the thermostats 
by their interaction with the system) due to the evolution 
is naturally defined in terms of the phase space contrac- 
tion: it is therefore interesting to see that even in absence 
of friction entropy production occurs and actually it can 
be identified, in the thermodynamic limit, with the same 
quantity that would arise in thermostats realized via ar- 
tificial frictional forces. The latter are widely studied in 
the numerical simulations as approximations to infinite 
systems in a thermodynamic limit, because it is not pos- 
sible to simulate really infinite systems. See Sec.V 

An important question is whether time evolution 
changes the configuration x into S t x but keeps the temper- 
atures and densities of the thermostats constant at least 
with fio -probability 1 and for any finite time. This is part 
of the more general question whether the spatial average 
of an intensive observable remains constant in time. 

A simple, partial but quantitative, formulation is in 
terms of the number Nj^iStx) of particles of S t x, of 
the kinetic energy Kj^(S t x) and of the potential energy 
Uj,h{Stx) of the configuration Stx into which x evolves 
at time t, inside a ball A centered at the origin. Consider, 
then, Vj > ,the limits (if existent) 

l im ( N ^(Stx) U jAn (Stx) K jiAn (S t x) \ 

n^ooV |A n n%| ' |A„n%| ' |A„nfij| /' 1 ' ' 

Under the above "no phase transition" assumption on fi 
we shall prove: 

Theorem 1: The limits in Eq.1.8 exist with ^-pro- 
bability 1 for all times and are time independent. The 
limits will be respectively Sj,Uj and ^SjksTj with fi - 
probability 1, as in Eq.1.3. 

Remark: This shows that the thermostats keep, in the 
thermodynamic limit, the same temperature and density 
that they had in the initial state: a property that has to 
be required for the model to adhere to the physical intu- 
ition behind the empirical notion of thermostats. Hence 
density and temperature of the thermostats are constants 



of motion. We shall show that more generally many other 
intensive observables are also constants of motion. 



II. INTENSIVE OBSERVABLES 

The definition of an hv -particles intensive observable 
is in terms of a smooth function T(q 1 ,q 1 , . . . ,qh,qh) on 
R 2dh vanishing for h ^ hr and which is "translation 
invariant", and with "short range" rr. 

This means that T = if the diameter of X = 
(qi, ...,qh) exceeds some ?t > and, denoting by 
T(:(X, X) the configuration (qx + £, q lt ...,q h + £, q h ), it 
is T{T t {X,X)) = T(X,X), V£ g R d . 

Given a region W the function Gw of x = (X, X) 

Gw{x) d = J2 ( 2J ) 

Ycxnw 

defines a "local observable" in W C R d with potential V. 

We shall say that Gw is an observable of potential type 
if T(Y,Y) depends only on Y, while if it depends only on 
Y it will be called of kinetic type. 

Then, if V„ d = f nj n A n , \V n \ d = f volumc(K), 

Definition 1: The "local average" of T on the configu- 
ration x = (X,X) is \V n \~ 1 Gv n (x). The corresponding 
"intensive observable" in the j-th thermostat is 

g{x) d g Um J-G Vn (x), (2.2) 

rwoo \ V n \ 

if the limit exists. Furthermore, given /i , define the "in- 
tensive fluctuation" of G (in the j-th thermostat) 

A G (z) lim ( -L GV>)- w (-LGyJ) (2.3) 

™^°° \v n \ \v n \ 

d = lim A G y n (a;), 

n— ¥oo 

if the limit exists. 

Remark: The notation requires keeping in mind that Gy n 
depends also on j (because V n — flj fl A„): however for 
simplicity of notation the labels j on V n and Gy n will not 
be marked. 

Properties of intensive observables can be derived from 
various assumptions on the initial distributions of the 
particles in the various regions f2j which, we recall, are 
distributed independently over j = 1, . . . , v and depend 
on the v pairs of parameters f3j , Xj . 

The simplest assumption is perhaps the uniqueness of 
the tangent plane to the graph of the pressure in var- 
ious directions, which could for instance be insured by 
the uniqueness of the translation invariant states of our 
particles system with parameters f3j , \j . 
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Let G be an observable of potential or kinetic type; 

def 

and suppose that iJ ,A,r( x ) = Hq,a(x) + 9G\(x) is su- 
perstable for \9\ small enough (i.e. there exist constants 
a > 0,6 > such that for all balls A it is H ^^(x) > 
aN 2 /\A\ - bN for all configurations x = (X, X) with N 
particles and with X C A and V \9\ < 9q for some 9 > 0. 
We call G an " allowed observable" . For such observables 
it is possible to define, for \9\ small, the "pressure" 



A-s-oo | V | 



Zj(0) 
Zj(0) 



(2.4) 



with given by Eq.1.2 with the energy 9Gy(x) 

added in the exponential. It is -P(O) = 0. 

It is important to stress that P(9) is, in the geome- 
tries in Fig.l considered here, independent of the special 
geometry considered for the f2 3 as long as the conical 
containers have d-dimensional shape (i.e. they contain 
balls of arbitrarily large radius) . 

In this context we can derive the following result: 

Theorem 2: Let G be an allowed observable of po- 
tential or kinetic type. If P(9) is differ entiable at 9 = 
0, then with [Iq -probability 1 the limit as \V n \ — > oc 
of j^—\Gv n (Stx) exists /j,Q-almost everywhere and is t- 
independent. 

Remarks: (1) The differentiability assumption of P(9) 
has the meaning of uniqueness of the tangent plane to 
the graph of the pressure p "in the direction of G" : such 
uniqueness is a "generic" property, see [7] for the lattice 
gas case. 

(2) The superstability of H -^(x) + 9G-^(x) is a very 
strong condition: it is certainly satisfied if 

(i) T(X,X) = 1 for \X\ = 1 and otherwise, or if 

(ii) F(X,X) = \q 2 for \X\ = 1 and otherwise, or if 

(iii) T(X, X) = unless X = (q,q') and in such case 
r(q,q')=<p(q-q'), 

therefore theorem 1 is a corollary of theorem 2. 

We also expect that the intensive observables will have 
very small probability of being appreciably different from 
their average values, and precisely a probability bounded 
above by an exponential of the volume |A„|. This will 
mean that the observable G satisfies a kind of large de- 
viations property: 

Theorem 3: Under the assumptions of theorem 2 the 
fio -probability that the fluctuation Ao,A n (Stx) differs 
from by more than e > tends to exponentially fast 
in \V n \ as n — » oc, Ve > 0. 

Remark: The assumptions in theorems 2,3 are satisfied 
by many observables in the Mayer expansion convergence 
region in the plane Xj , (3j , [8] . They are also believed to 
be satisfied quite generally for observables generated by 
a potential T. In particular they hold generically if T is a 



linear combination of the potentials (i),(ii),(iii) in remark 
(2) above. 

The proof of theorems 2,3 arc presented in Sec. IV. 



III. TIME EVOLUTION 

A quantitative existence theorem of the dynamics can 
be conveniently formulated in terms of the quantities 

to, r v and W,Af,v\, \\x\\\ defined as 



W(x;£,R) 



def 1 



E 



,mq~ 



ip(o) V / 2 



+\ Yl viQi - v) + ^(u) + vi )) i 



A/j(x) d = number particles within r v of £ G M. d , (3.1) 



def 

-x'iW = \qi - qtl/v! + \q t - q[\jr v 



Let log + z d = max{l,log 2 |z|}, g c (z) = (log + z) c and 

W(x-tR) 



cq(x) — sup sup 

6 R>gcd/r v ) 



R d 



Call H( the configurations in H with 



(3.2) 



(3.3) 



(2) 



(1) £ c (x) < oo 
Nil, A n ) U(j,A n ) K(j,A n ) 

|A„n%|' |A„n^-|' \A n nQj\ n ^°°' ^ ,U,J ' 2^ 



d5j 



with A„ the ball centered at the origin and of ra- 
dius 2 n r ip 5j,Uj,Tj, given by Eq.(el.8) if N(j,A n ), 
U(j, A n ), K(j, A n ) denote the number of particles and 
their internal potential or kinetic energy in Clj PlA„. Each 
set %c, has /io-probability 1 for £ > 1/d, [2, 9-11]. Then: 

Theorem 4: Let d < 3, then Hi/d has fio -probability 1 
and S t x exists for ^-almost all x <G 'H\/d an d Vt > 0. 

Given (arbitrarily) a time 6 > 0, if £ d = £\/d( x ), and 
I ft (0)| < 2 k r v there are c = c(£,@) < oo,c' = c'(£,Q) > 
such that Vn > k and \/t < 9 



\qi(t)\ < evi fca, 

distance(gi(i), 9(Ujf2j)) > c'r v k~i 
Af 5 {S t x) < ck 1 / 2 

IKStaOi-^aOiH < e- c ' 2 " /2 , n > k. 



(3.4) 



This is proved in [5, theorem 7] for d = 2 and in [6] for 
d = 3 (the latter reference covers also the case d — 2 via 
a somewhat different approach). 



■5 



Remark that the theorem does not state that the sec- 
ond of Eq.3.3 holds: in [5, 6] it is however proved, in 
addition to theorem 4, the weaker statement that the 
liminf of —fefpiTT^ is not smaller than 5 of the corre- 
sponding r.h.s.; and the same is true for the other two 
quantities in Eq.3.3. 

A corollary of the main results of this paper will be 
that the limit relations in Eq.3.3 will hold for all t > 0. 



IV. CONSTANTS OF MOTION 

Let T be an fr-points local observable of potential type, 
V n = Slj n A„. Under the assumptions of theorem 2 we 



Oj n A„. 

first show that lim n 

, def 



\Vn\ 
r log 



1 1™ / p -8G Vn \ 



- g exists. 

this is smooth 
is 

= 0, 



Define P n {6) - 
and convex in 9 and its unique derivative at 6 = 

9n = fp^yMo(GV n ); therefore, remarking that -P(O) 
it satisfies P n {0) > 9 g n . 

The limit P(9) as n — > oc of P n (6) is the same that 
would be obtained if V n was replaced by the full ball A„ 
and filled with particles at temperature (3J 1 and chemical 
potential Xj. 

Any convergent subsequence g ni defines therefore a co- 
efficient g with the property P(0) > 6g. Hence, by the 
assumed uniqueness of the tangent to P(9) at 6 = 0, it 
follows that g is uniquely determined thus implying that 

the limit g = lim^oo g n exists. 

Let g n = { \V n \~ 1 G Arl ) M and, given 7 > 0, let X E ,-y, n 
to be the set of points in Hi/d with £(x) < E, G An (x) < 
(g n + ^)\V n \ and which, under the evolution, reach in 
a time r 7i „(x) < 6 and for the first time, a point of the 
surface 



def 



{x\\V n \- 1 G An (x) = (g n + 1 )}. (4.1) 



If for all E and for all small 7 > it is J2 n u o {^E,j,n) < 
+00 then it will be limsup^^^ |T^ l | _1 GA„(<5 , *a;) < g, 
with u -probability 1 (by Borel-Cantelli's estimate); 
changing T into — V it will follow, again with ^ 
probability 1, that the liminf is > g: notice that the 
change in sign of T is possible by the conditon on G to 
be an "allowed observable" , as introduced before Eq.2.4. 

This remains true if for all small 7 there exists j n £ 
[7,27] such that J2 n ^o(^E nrl ,n) < +00. 

If x £ XE,-y,n the phase space contraction, when phase 
space volume is measured by u , within time t is, [5, 6], 



s(x, t) 



3>0 



(4.2) 



Fj,L (t) = / Xo .#(X (t)). 



where Q j (t) d = Xj (t) 

By theorem 4, L (t) is uniformly bounded as n — > 00, 
for < t < 9, by the first of Eq.3.2, by a quantity C 
(only depending on E, no, 0). 



Therefore by a quasi-invariance lemma, [2, 12], [5, Ap- 
pendix H], the probability ^oi^E.-y+e^) can be bounded 

Ve e [7, 2 7 ] by 



\G\ , r G An (x) 



\V n 



(gn+J + s)) (4.3) 



where G denotes the time derivative (at t = 0) of 
G An (S t x) (to be computed via the equations of motion) 
evaluated on the surface E„ i7 + e , see Eq.4.1. 

Integrating Eq.4.3 over efe/7, fi (^E,n,i n ) can be 
bounded by 



^ / Mo(^)M x(7 < °jM ^ gn < 2 7 ), (4.4) 

7 J \V n \ \V„\ 

with G = J2xcv n 'Eqex d q T ( x ) 9- B y Schwartz' in- 
equality 



G2I- 1 uo({x : 7 < %^ -9n< 2 7 }) 1/2 (4.5) 



\V n \ 



because from Eq.2.1 for T 



/i (G 2 ) 1/2 < Ci\V n \ 



(4.6) 



obtained via superstability bounds, using the Maxwellian 
distribution for q. 

The probability in Eq.4.5 is bounded above by Chebis- 
chev inequalities (quadratic or exponential) by both av- 
erages 

J <*±f , (GAAx)/\V n \-9n) 2 \ 

Ie d g {e 0(G An -\V n \ {gn+1 ))^ o (47) 

V9 > 0. This implies the existence of 7„ e [7, 27] with: 

Mo(^,„, 7 J < C 37 - 1 J(n), J(n) 2 = I,I e (4.8) 

Therefore we look for assumptions on the thermostats 
structure (i.e. on \j,/3j,(p) under which J(n) tends to 
zero fast enough making ^ n Ho{^E,n,j n ) < 00. In this 
case theorem 2 will follow from Borel-Cantelli's lemma 
and the arbtrarincss of 7. 

As a consequence of the above bounds, basically fol- 
lowing from the uniqueness of the tangent plane in the 
direction T, the proof of theorem 2 can be completed as 
follows. Fix 7 > and remark that 

Ie = (e eUr - v " ) e - e (sn+7)|v„| < e -6~/\V n \+ri(9,V n ) ( 4 g) 

Continuing the argument leading to the existence of the 
limit of g n , at the beginning of the section, the correction 
term rj(9,V n ) is bounded, cas follows: 
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(a) ^log( e ^,v„ is P n {9) - p„(0) (notice: P n (0) 

= 0) and converges to P{9) — P(0) as V n — > oo for \9\ < 
6q, if 9q is small enough so that the potential ip + fij 1 9 
is superstable V \9\ < 0$, j = 1, . . . , v. By superstability 
the limit exists for \9\ < 9 and it is a limit of functions 
P n {9) which are convex for \0\ < 9 . Hence the limit is 
uniform: \P{9) - P n {9)\ < o(\V n \) for |0| < 9 , 

(b) the g n in the exponent in Eq.4.7 has just been shown 
to be g n \V n \ = g\V n \ + o(\V n \), so that -0 g n converges 
to —9 g with an error 9 o(\V n \), 

(c) (P{9) - P(0) -9g)\V n \ is (by the uniqueness of the 
tangent plane) o{9) \ V n \. Hence 



V (9, V n ) - 7 0„|K| < -\l9 n \V n \ (4.10) 

+ ( - + OJ wf + °M ^ -h 6 n\V n \ 

and choosing 9 n tending to so slowly that the exponent 
of the r.h.s. of 4.9 tends rapidly to oo, for instance if 
9 n = max( TS |^, ^^T^fp), we see that I 0n -^^0 so fast 
that Ho(^E,n.'y n ) is summable in n implying theorem 2 
and of its special case theorem 1. 

Theorem 3 also follows from the existence of the limit 
for g n because Ig yields a summable bound on J, hence 
on p Q (A 2 GA J. 

Remarks: (1) Uniqueness of the tangent plane can be re- 
placed by assumptions on the decays of correlations in 
the distribution p somewhat stronger than just requir- 
ing its extremality among the DLR distributions in the 
geometry in Fig.l. 

(2) Sufficient estimates can be formulated as follows: 
Pj(xi, . . . , x n ) be the n-points correlation function in 
the j-th container: by superstability pj < C n , [9]. If 
x = (q,q) and £ G flj, extremality of po, implies, [4, 9], 
for x\ , . . . , x n and yi, . . . , y m ) with positions in Qj : 



r. We could actually consider more general observables 
of the form (in Vi j ) 



\Pj(xi,- ■ ■,x n ,T^y 1 , . . .,Ttfj m ) 

~Pj{xi, ■ ■ . ,x„)p J -(r e y 1 , . . . ,T^y m )\ 



{— s-oo 



.11) 



Assume that Eq.4-.11 holds in the stronger sense that 
the l.h.s. is bounded by r]R^ m ,n{£,) if the positions of 
X\ , . . . , x n and yi, .. . , y m ) can be enclosed in a ball of 
radius R. 

Theorem 5: // there is a constant Cn tTn , n < oo such 
thatVR,m, n (0 < C R , m ,n\Z\- a{R ' m ' n) with a(R,m,n) > 
and if liixiA->oo t^-iP<o{Gv„) = 9 exists, then limA^oo 
A G y n (x) = and lim A ^oo \VZ\ G v n (S t x) = g with p 
probability 1. 

Remarks: (1) Thus if po has a power law cluster property 
all intensive observables admitting an average value, over 
space translations, at time are constants of motion. 
(2) With the above assumptions we avoid use of the expo- 
nential Chebishev inequality and we may thus drop the 
superstability condition in the definition of the potential 



lim 



1 



|% n A„| J re m. d :T r Acn : ,nA„ 



T r f(x)dr (4.12) 



where / is a cylindrical function in A (i.e. it does not 
depend on the particles outside A) and r r denotes trans- 
lation by r. If the power law cluster property is satisfied 
and po a.s. the limit in 4.12 exists at time 0, then the 
intensive observables 4.12 are constant of motion under 
the assumption that / is smooth and grows at most poly- 
nomially with the number of particles. 
(3) The assumption certainly holds in the cluster expan- 
sion convergence region, [3] and [13, Sec. 5. 9], i.e. high 
temperature and low density, without extra assumptions. 

Proof: Consider the first of Eq.4.7 and choose 7 = j n — 
i. The numerator tends to as \V n \-< R > n ' n )/ d if the 
potential T for the observable Gy n vanishes when the 
diameter of the set {x\, . . . , x n } exceeds R. 

The estimate Eq.4.8 implies that j^—^Ag,v„ (Stx) tends 
to with /xo-probability 1 for all t < kt with k integer 
and t n > (arbitrarily fixed). Hence if the average of 
roMo(GV n ) exists it exists for all times and has a time- 
independent value. 

V. ENTROPY AND THERMOSTATS 

Entropy production rate (due to the action of the sys- 
tem upon the thermostats and identified with the rate 
of their entropy increase, which is finite even though 
the thermostats entropy is infinite because the ther- 
mostats are infinite) is defined in terms of Qj = — ~Kj ■ 
dxjUo,j(X-c Xj), which is the work per unit time, per- 
formed by the test system on the j-th thermostat. Since 
Qj is interpreted as the heat ceded by the system to the 
thermostats the entropy production in the configuration 
x is given by er (x) = Ej>o PjQA x )- 
If the volumes in phase space are measured by the distri- 
bution po this quantity differs from the contraction rate 
of the phase space volume by /3o(Qo + £o) = flo(Ko + Uo) 
and K + Uo is expected to stay finite uniformly in time. 
If so the statistics of the long time averages of the phase 
space contraction rate and of the entropy production rate 
will coincide (however this is not proved as the theorems 
above only concern what happens in a arbitrarily prefixed 
but finite time interval) . 

In other words in the frictionless thermostats model and 
in the isoenergetic thermostat models, [6], the entropy 
production can be identified with the phase space con- 
traction, possibly up to a time derivative of a quantity 
expected to be uniformly finite in time. Furthermore the 
entropy production is the same in both models of ther- 
mostats if the thermodynamic parameters of the ther- 
mostats (Sj,Tj, j > 0) are the same: this follows from 
the equivalence theorem between frictionless and isoen- 
ergetic thermostats, [6, theorem 1], which states that un- 
der such conditions the microscopic motions of the two 



7 



models starting from the same initial condition remain 
identical forever with /i -probability f . 

Other thermostats can be considered: for instance the 
isokinetic thermostats. At a heuristic level analogous 
conclusions can be reached, [14]. 

Considering external thermostats as correctly repre- 
senting the physics of the interaction of a system in con- 
tact with external reservoirs has been introduced in [151. 



Their analysis was founded on the grounds of 

(1) identity, in the thermodynamic limit, of the evolution 
with and without thermostats 

(2) identity of the phase space contraction of the ther- 
mostatted systems with the physical entropy production 
(up to a time derivative). 

For a more mathematical view see [14]. 
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